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I. Introduction 


An intense solar microwave millisecond spike emission (SMMSE) event has been observed on May 16, 1981 
by Zhao and Jin at Beijing Observatory (a> = 2 tt x 2.84 GHz)W. The peak flux density of the spikes is high 
to 5 x 10 5 s.f.u., and the corresponding brightness temperature (BT) reaches - 10 15 K. 

In order to explain the observed properties of SMMSE, we propose in this paper that a beam of electrons 
with energy of tens kev injected from the acceleration region downwards into an emerging magnetic arch 
forms so-called hollow beam distribution and causes electron-cyclotron maser (ECM) instability. The growth 
rate of second harmonic X-mode is calculated and its change with time is deduced. It is shown that the satura- 
tion time of ECM is t s * 0.42 ms and only at last short stage (At < 0.2 t s ) the growth rate decreases to zero 
rather rapidly. So a SMMSE with very high BT (T^ > 10 15 K) will be produced if the ratio of number densi- 
ty of nonthermal electrons to that of background electrons, n s /n e , is larger than 4 x 10 -5 . 

II. Model and Excitation of X-mode 


If the induced electric field is larger than the 
Dreicer field, it will be able to accelerate electrons 
along the field to energy of tens kev in one milli- 
second. The injected electrons with a large pitch 
angle gyrate around the magnetic field B and form 
so- called hollow beam distribution (Fig. 1) as fol- 
lowing equation^ 21 : 


Fty.v,,) = (tt 3/2 dj d„ v 2 ^,,)- 


exp 


|-(*3 - (S31 



Figure 1 . Model of hollow beam driven ECM 
for SMMSE. 
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It drives ECM and we have the growth rate at the m-th harmonic 121 
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c is light velocity, n is refraction index satisfying the Appleton-Hartree dispersion relation, a is the angle be- 
tween wave vector K and B, « is wave frequency, 
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8. 


Assuming that the energy of nonthermal electrons E = 60 keV, pitch angle <t> - 60°, i.e., v 0i /c = 0.386, 
v 01! /c = 0.223, and dj^ = 0.2, d n = 0.4, « p /8 e = 0.4, it is found that the second harmonic (m = 2) X- 
mode will grow fastest if to = 2.02 8 e , a = 65°, and the growth rate is 


a,® = 0.096— 8. 

n. ‘ 
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m. Saturation of ECM 


The excitation of ECM leads to diffusion of energetic electrons in velocity space, which in turn weakens the 
instability. According to quasilinear theory we have 131 


i^V’Vn) - ^ 4; F(V 1 ’ V "> 

W,v L ,v„) =^r ydKSiT(t) 6 (« - ^ - K„ v„) 


(3) 

(4) 


-|g?(t) - 2« ( (K,t)fit(t) (5) 

where 7 = (1 -(v 2 + v n 2 )/c 2 ) - *, D(t, ,v n ) is the diffusion coefficient and S-£the wave spectral density. 
Let 


r 


=f dt'D(t', v i , v„) 


( 6 ) 


we can rewrite Eq. (3) as 

3F1 d j>F 

dr _ v £ d Vj^ i J » 1 

then obtain the distribution for t > 0 (i.e. t > 0). 

F(r, Vi ,v„) dy F(0,y, V|1 ) e~ h(^)y 

where F(0, , v n ) is the initial distribution of Eq. (1), I Q (x) is modified Bessel function of zero order with argu- 

ment x. 



Figure 2. Plots of <o< 2) /fl e •rig/ng: (a) vs r, (b) vs time t, all for X-mode with parameters: E = 60 keV, 
0 =60°, a = 65°, ojp = 0.412 e , d^ = 0.2, d u = 0.4 and o) — 2.02Q e . 
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Substituting Eq. (7) into Eq. (2) and using the following relations 

(JF \ _ /_3F \ /dF\ (Jt\ 

\ av,./ y av J' yaT /’ 11 ya vj' 

/ dF \ r /3F\ 

*V dv n/ r V n W V| ‘ 

and Eqs. (4)-(6), we obtain the change of u\ 2) with r (Fig. 2a) and then with time t (Fig. 2b), which gives the 
saturation time of ECM t s (r = r s ) corresponding to o)\ 2) — 0. We have t s *= 0.42 ms and find out that wf 2) 
almost doesn’t change until t * 0.8 t s and only at last short stage (At < 0.2 t s ) it decreases to zero rapidly. 

In addition we have obtained the following relations 

t - J lJ 2 m\ 

s 2« io V W J (8) 

~4^ TsWioK '' V °" (9) 


where W(t s ) and W 0 are the saturation and initial energy density of the wave, respectively, u>j 0 is the initial 
growth rate and K* = rjulc cosa. 

The above calculation is rather complicated. For simplicity we give an evaluation for saturation value r s . 

Taking only the part of distribution (1) containing v^ and introducing x s v j/c, v = v 0 /c, 

A — d|V, we have 


F 0 (x) = (^rc 2 A 2 )- 1 exp 


i - mi 


a^ = _ 2(x-v) 
d\ ttc 2 A 4 


^ exp rrr 


X-vV 


If X = Xo = _ dj_/V2) v, then 3F 0 /dx > 0 (so wj > 0) and has its maximum. Introducing T = t/c 2 , 

X = Xq’ z =y/c, we have for t > 0 (t > 0) 
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we can divide the integration in Eq. (10) into two parts: 
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where z 0 ^ lOT/x and in Q, (0 < z < z 0 ) the Bessel functions I 0 , I, have been expressed approximately for 
small argument but in Q 2 (z G < z < 00 ) for large argument. It is easy to see that in general we have 
Q, < 0 and Q, > 0. 

If z 0 « x? i e. T « 0. 1 x 2 , Q2 is dominant and d¥ld\ > 0. This is the case at the start of the growing wave. 

If z 0 >> x> i e. T >> 0.1 x 2 > Q, becomes dominant, so dF/d\ < 0. This is the case after saturation. 

If z 0 - X* i-e., T - 0.1 x 2 , we have Q 2 > | Q, | and dFIdx > 0, but | Q, | gets larger. This corresponds 
to the case that the growth rate decreases obviously. 

Hence we can take T s =0.1 v 2 or r s =0.1 v 2 ^ as the saturation value when Q, can be comparable with 
Q r For the chosen parameter, v 0 /c = 0.386, we have T s = 0.015 which is in agreement with numerical 
calculation (T s * 0.013). 

IV. Gyroresonance Absorption and BT of Radiation 


The gyroresonance absorption (GA) is a serious problem for maser radiation. The optical depth of the mth 
layer (m > 2) for x-mode is HI 


coV rn 2m >/2 co L b / T e W 1 
2 m_l rn! c ym, c 2 J 


with L b =B dB 
dL 


For co= 27 t x 2.84 GHz and m =3 we have 




K 

10 8 


Because L B > 1000 km in general, an intense maser radiation (T^ > 10 l4 K) can hardly be received unless 
T e < 2xl0 6 K and cop/co < 0.1 

It seems possible, however, that some magnetic arches with smaller dimension and higher magnitude may 
emerge from an active region where the magnetic configuration is more complicated, then it may be true that 
L b < 1000 km (e.g. L B = 500 km) at larger angle to B and the third GA will be weaker, so the radiation 
may escape from the corona. 
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The BT (T 5 ) and flux density (S) received are 



where W is the energy density of the wave, R = 1.5 x 10 13 cm is the distance between the Sun and the 
Earth, A is the surface area of the radiation source, which is about 10 15 cm 2 for SMMSE, Aa> is the bandwidth 
of the radiation (Aoj/co « 0.01 according to Droge 151 and our calculation), AQ r is the radiation solid angle. It 
is easy to see from Fig. 3 that the propagation angle a will be from 60° to 70° if the pitch angle <t> = 60°, 
hence Aft r = 1. 



Figure 3. A plot of o>pVfl e • n e /n s vs <j> for Figure 4. A plot of v^ part of distribution of nonther- 

different a: 55°, 60°, 65°, 70° and 75°. mal electrons for different T with x = v^/c, T 

Other parameters are the same as in Fig. 2. = r/c 2 , v = v 0 ^/c = 0.386, and x 0 = 0.86v. 


Assuming n = 0.97, a = 65° (for the fastest growing second harmonic of X-model) and n s /n e = 
4 x 10" 5 , r s /c 2 = 0.013, W(t s )/W 0 = 2 X 10 10 , e ~ r 3 = 0.1, we have obtained from Eqs. (8)— (13) 

t s - 0.36 ms, W(t s ) * 5.9 x 10 -5 ergs cm -3 
T b * 5 x 10 15 K, S ~ 5.5 x 10 5 s.f.u. 

These results are in agreement with observations. 
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V. Summary and Concluding Remarks 

It has been shown that ECM is an effective mechanism to produce electromagnetic radiation. Its BT may be 
higher than 10 15 K and it can escape from the solar corona without any transformation. The simple procedure, 
high efficiency, short time scale, strong directivity of radiation, and high degree of circular polarization are its 
advantages over the plasma radiation mechanism to explain SMMSE. It is calculated in this paper that 
SMMSE with very high BT (T^ * 5 X 10 15 K) would be produced by hollow beam driven ECM if E = 

60 keV, 4> = 60°, a = 65°, n s /n e = 4 X 10" 5 , w p /Q e = 0.4, L B = 500 km, and T e = 2 X lO^C. 

It is shown that the range of pitch angle <t> for exciting ECM is very wide (> 30°, see Fig. 3) and the 
energy of nonthermal electrons may be from ~ 10 keV to 100 keV or higher. Besides the hollow beam distri- 
bution, other velocity anisotropy as a loss-cone distribution can also excite ECM under appropriate conditions. 
Hence it seems that ECM is often produced above solar active regions. The main reasons why there appear 
rarely intense SMMSE events may be due to the powerful directivity of radiation and strong third GA. 

It is easy to see from Eq. (7) and Fig. 4 that with the evolution of the energetic electron’s distribution, the 
transverse velocity component of them is decreased but the parallel one isn’t. As a result, the pitch angle is 
getting smaller, so the electrons, after losing most of their transverse energy, get through the magnetic mirror 
points and rush into the transition region to produce hard X-ray bursts, so with which the SMMSE correlate 
very well 161 . 
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